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Abstract 



We derive the deformed s^(2) Gaudin model with integrable boundaries. Start- 
Q\ | ing from the Jordanian deformation of the SL(2)-invariantYangR-matrix and generic 

solutions of the associated reflection equation and the dual reflection equation, the 
T^j- \ corresponding inhomogeneous spin-V2 XXX chain is obtained. The quasi-classical 

expansion of the transfer matrix yields the deformed s£(2) Gaudin Hamiltonians 

with boundary terms. 
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I Introduction 

Gaudin models have applications in many areas of modern physics, from quantum 
optics [HEI to physics of metallic nano-grains, see |3| and reference therein. A model of 
interacting spins in a chain was first considered by Gaudin HH5]]. In his approach, these 
models were introduced as a quasi-classical limit of the integrable quantum chains. 
Moreover, the Gaudin models were extended to any simple Lie algebra, with arbitrary 
irreducible representation at each site of the chain |5ll . 

The rational si (2) invariant model was studied in the framework of the quantum 
inverse scattering method |6|. In his studies, Sklyanin used the si (2) invariant classical 
r-matrix. A generalization of these results to all cases when skew-symmetric r-matrix 
satisfies the classical Yang-Baxter equation ]7] was relatively straightforward |8j. There- 
fore, considerable attention has been devoted to Gaudin models corresponding to the 
the classical r-matrices of simple Lie algebras |9] and Lie superalgebras IfTOi ni. 12j. In 
the case of the s£(2) Gaudin system, its relation to Knizhnik-Zamolodchikov equation 
of conformal field theory HT3l [Ml [15l or the method of Gauss factorization [16], pro- 
vided alternative approaches to computation of correlation functions. The non-unitary 
r-matrices and the corresponding Gaudin models have been studied recently, see |T7fl 
and the references therein. 

The quantum inverse scattering method (QISM) ffl8l [19] |20l as an approach to con- 
struct and solve quantum integrable systems has lead to the theory of quantum groups 
|2Tll22| . A particularly interesting feature of quantum groups is a transformation that 
is called twist [23| and allows constructing new quantum groups from already known 
ones. Although the twist transformations generate an equivalence relation between 
quantum groups, they produce different R-matrices. These new R-matrices can in turn 
lead to new integrable systems |24| . 

Twist of a quantum group, or more general Hopf algebra A, is a similarity transfor- 
mation of the coproduct A : A — > A <g> A by an invertible twist element 

r(l) « r(2) 



T = Y^f- } ®fj eA®A, 



A(fl) -> A t (fl) = TA(a)T~ l , a e A. (1.1) 

In order to guarantee the coassociativity property of the coproduct, the element T has 
to satisfy certain compatibility condition, the so-called twist equation [23 1 



where 



Tn (A ® id) (F) = J- 23 (id ® A) (T), (1.2) 
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The transformation law of the coproduct also determines how the corresponding uni- 
versal R-matrix is changed 



K -> ftW = F 2 iKT-\ T 1X = J2fj 2) ® fj 1] - (1-3) 



/2)^(l) 

This new R-matrix allows us to build and study new integrable models [24 J. 

A particular solution of the twist equation is provided by the Jordanian twist ele- 
ment for the universal enveloping algebra U(s£(2)) of the s£(2) Lie algebra [251I26H . It 
was extended to the s£(n) case in |27ll28ll . The Yang R-matrix determines the Yangian 
y(s£(n)), see HQ. The Heisenberg XXX spin chain is related to the Yangian y{s£{2)) 
and the universal enveloping algebra of s£(2) is a Hopf subalgebra of the Yangian, 
U(s£(2)) C y(s£(2)). It can be shown that the transformation of the Yang R-matrix 
by the Jordanian twist element yields the R-matrix of the twisted Yangian yg(s£(2)) 
||30l[3Tl . Because the twist preserves the regularity of the R-matrix the Hamiltonian of 
the deformed Heisenberg XXX spin chain with periodic boundary conditions can be 
calculated I30|. Although it can be seen that the extra terms added to the XXX Hamil- 
tonian do not change the spectrum of the model, the explicit form of the Bethe states is 
not obvious. 

The deformation by the Jordanian twist of the XXX model with non-periodic bound- 
ary conditions was studied in 11321 . A way to introduce non-periodic boundary condi- 
tions compatible with the integrability of the quantum systems solvable by the QISM 
was developed in [33|. The boundary conditions at the left and right end of the sys- 
tem are expressed in the left and right reflection matrices. The compatibility condition 
between the bulk and the boundary of the system takes the form of the so-called re- 
flection equation |34ll35l . The compatibility at the right end of the model is expressed 
by the dual reflection equation. The matrix form of the exchange relations between 
the entries of the Sklyanin monodromy matrix is analogous to the reflection equation. 
Together with the dual reflection equation they yield the commutativity of the open 
transfer matrix |33ll36ll . The general solution of the reflection equation associated with 
the Jordanian deformation of the SL(2)-invariant Yang R-matrix was given in t32l , as 
well as the Hamiltonian of the deformed XXX spin chain with the general boundary 
terms. 

Here, we study the deformation by the Jordanian twist of the s£(2) Gaudin model 
with non-periodic boundary conditions. The starting point to obtain the Gaudin model 
in the framework of the QISM is the monodromy matrix of the corresponding inho- 
mogeneous spin chain [6|. The quasi-classical expansion of the transfer matrix of the 
periodic chain, calculated at special values of the spectral parameter, produces the 
Gaudin Hamiltonians 11371 l38l . This approach was generalized to the case of non- 
periodic boundary conditions in ||39l . These results were later extended to non-diagonal 
reflection matrices |40[ and to other simple Lie algebras |4T1 . Essential steps in the 
derivation of the Gaudin Hamiltonians with boundary terms are some normalization 
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conditions which have to be imposed both on the R-matrix and on the left and right 
reflection matrices. In general, these normalization conditions are not relevant in the 
study of the spin chain. However, they enable the quasi-classical expansion of the trans- 
fer matrix which yields the corresponding Gaudin Hamiltonians with boundary terms 

nana. 

This paper is organized as follows. In Section II, the Jordanian deformation of the 
SL(2)-invariant Yang R-matrix is reviewed. The general solutions of the correspond- 
ing reflection equation and the dual reflection equation are given. In Section III, the 
Jordanian deformation of the inhomogeneous XXX spin chain with N sites is studied. 
It is shown how the quasi-classical expansion of the transfer matrix yields the Gaudin 
Hamiltonians. These Hamiltonians are calculated in Section IV. 

II Twisted Yang R-matrix and reflection equation 

We start this section by reviewing the Jordanian twist element. The three generators of 
the s£(2) Lie algebra are h, X ± , such that 

[h,X ± ]=±2X ± , [X + ,X-] = h. (II.l) 

The universal enveloping algebra U(s£(2)) admits the Jordanian twist element ||25ll26ll 



7 = exp(h <g> ln(l + 29X+)) = e hm G U(si(2)) ® U(s£(2)) (11.2) 

which satisfies the following equations 

(A ® id) (T) = F13T23, (0.3) 

(id ®A t )(T) =J r i 2 Ji 3 . (0.4) 

The coproduct A is the usual coproduct of U(s£(2)) and A t is the twisted one. In order 
to check the equations above, it is important to notice that h is primitive with respect to 
A and that a is primitive with respect to A t , 

A(/i) =h®l + l®h, (0.5) 

A t (cr) =crtg)l + l(g)cr. (II.6) 

The equations 011.31 ) and (|II.4|) imply the twist equation ((12)) , 

Tn (A ® id) (T) = TnTnT-x, (0.7) 

J23 (id ® A) (T) = J23 (id ® A) ( T)T^ ■ J23 = (id ® A t ) (T) J" 23 = ^12^13^23- 

(0.8) 



Therefore the Jordanian twist element ( |II.2| ) satisfies the Drinfeld twist equation (IL2" 
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The Heisenberg XXX spin chain is related to the Yangian y(s£(2)) and the SL(2)- 
invariant Yang R-matrix 



R(A) 



/ 1 








A 




'7 



A + rj A + t] 
t] A 



0\ 







A + tj A + rj 
V 1 / 



(II.9) 



where A is a spectral parameter, rj is a quasi-classical parameter. 

The matrix form of the Jordanian twist element pi.2|) in the spin Vi representation p 

is F12 E End (C 2 <g> C 2 ) 



F 12 = (p ® p) jr = exp (a z <8) 0cr+) = 1 + 0cj z <g> tr 4 



/ 1 \ 

10 

1-0 

V 1 ) 



(11.10) 



where c" z , a ± = (o~ x ± ic 37 ) /2 are the Pauli sigma matrices. Hence, the twisted R-matrix 
has the form 



/ 



R^(A) = F 21 R 12 (A)F- 1 



-A0 



A 6 



A9 1 \ 



A + T] 

A 


A + T] 


A + T] 

-A6 


A + rj 


A + T] 

A 


A + rj 
A6 







A + rj A + rj A + rj 

\ 1 



(11.11) 



/ 



This R-matrix is also a solution of the Yang-Baxter equation 

R 12 (A - f /)R 13 (A)R 23 (^) = R 23 (f/)R 13 (A)R 12 (A - p), 



(11.12) 



we use the standard notation of the QISM lfT8l[T9ll20l to denote spaces Vj,j = 1,2,3 on 
which corresponding R-matrices R,y, ij = 12, 13, 23 act non-trivially. In the present case 
V\ = V2 = V3 = C 2 . In what folio ws we will only use the twisted R-matrix (|II.11|) and 
in order to simplify the notation we will drop the symbol (/) in the superscript. 

The Gaudin models are related to the classical r-matrices |6|. Therefore it is essential 
that, after setting 6 = —^rj, the R-matrix has the quasi-classical property B3l 



R(A,^,0)| 9= _^ = 1 + ^(A) + O^ 2 ), 



(11.13) 
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here r(A) is the corresponding classical r-matrix 



r(A) 



-11. 
x x f 

1 -1 

A X -« 
V / 



(11.14) 



which has the unitarity property 

J*2i(-A) = -ri 2 (A), 
and satisfies the classical Yang-Baxter equation 

MA),r 23 (^)] + [r 12 (A - f/),r 13 (A) + r 23 (f/)] = 0. 



(11.15) 



(11.16) 



Moreover, for the purpose of deriving the Gaudin Hamiltonians, it is necessary that the 
R-matrix (|HTT1) is normalized so that 13711551 



R(0,t])=V, 



(11.17) 



where V is the permutation matrix in C 2 <g> C 2 . The R-matrix (|II.11)> has the unitarity 
property 

R 12 (A)R 21 (-A) = 1, (11.18) 

but the PT symmetry is broken 



R 2 i(A) ^ Ri 2 W hh , 



(11.19) 



where R 2 i(A) = 7 3 Ri 2 (A)'P / and the indices t\ and f 2 denote the respective transposi- 
tions in the first and second space of the tensor product C 2 ® C 2 . The R-matrix also 
does not have the crossing symmetry but it dos satisfy the weaker condition 



{{{R 12 (A) t2 }- 1 } f2 }- 1 = #(A)M 2 R 12 (A + 2t])M 2 \ 



(11.20) 



. , .,. A(A + 2n) 
with g(A) — —± — : — ^r and the matrix 



(\+ v y 



M 



We note that a more general matrix 



M 



1 -26 
1 



1 a 
1 



(11.21) 



(11.22) 
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commutes with the R-matrix, 



M<g>M,R(A) 



0. 



(11.23) 



A way to introduce non-periodic boundary conditions compatible with the integra- 
bility of the base model was developed in [33|. The boundary conditions at the left 
and right sites of the system are expressed in the left and right reflection matrices K~ 
and K + . The compatibility condition between the bulk and the boundary of the system 
takes the form of the so-called reflection equation f34ll35| . For the left reflection matrix, 
it is written in the form 

R 12 (A - ^)X-(A)R 21 (A + y)K 2 (y) = K 2 (y)R 12 (A + ^K^(A)R 21 (A - y). (11.24) 

The compatibility at the right site of the model is expressed by the dual reflection equa- 
tion 



A 12 (A - y)K+ t (A)B 12 (A + y)K+\y) = K+ l (y)C 12 (A + y)K+ t (A)D 12 (A - y). (11.25) 

where the matrices A, B, C, D are obtained from the R-matrix of reflection equation 
(UnUlas 



Aiz(A) = (Ri 2 (A) f «) 1 = D 21 (A), 

B 12 (\)=((R t 2 \(A)Y 1 ) 2 = C 21 (A). 



(11.26) 
(11.27) 



Using property (III.20D , we can write dual reflection equation dll.251) in the equivalent 
form 



R 12 (-A + y)K+(A)M 2 R 21 (-A - y - 2^)M 2 l K+{y) ■- 
K+{y)M 1 R 12 {-A -y- 2^)M 1 - 1 X+(A)R 21 (-A + y). 

It can then be verified that the mapping 

K+(A) =K.-{-A-rj) M 



(11.28) 



(11.29) 



is a bijection between solutions of the reflection equation and the dual reflection equa- 
tion. 

The general solution of the reflection equation associated with the Jordanian defor- 
mation of the Yang R-matrix is 11321 



/ 



K-(A) 



d(A) 



C-A 



-A 2 



<pA 



tpA 



n 



) 



(11.30) 
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with arbitrary parameters £, <p and t/> and the function 



d(A) = £ - —A 2 + A^l + </>!/? . (11.31) 

The left reflection matrix K~ is normalized so that 

jr(A)K~(-A) = 1. (11.32) 

Finally, by setting 9 = —^rj we achive that matrix R~ (A) does not depend on the quasi- 
classical parameter rj, 

^1=0. (11.33) 

The right reflection matrix R + (A) is obtained by substituting (III30I) into dll.291) while 
keeping 6 = — £,rj. It is important to notice that 

limX + (A) ) K~(\) = 1. (11.34) 

7-X) J 

In general, the normalization conditions dll.321) and (III.34D are not essential in the study 
of the open spin chain. However, together with ( 1II.13[) and flII.17[) they enable the quasi- 
classical expansion of the transfer matrix which yields the open Gaudin model. 

Ill Deformed inhomogeneous XXX chain 

In this section we study the Jordanian deformation of the inhomogeneous XXX spin 
chain with N sites, characterized by the local space Vj = C 2 and inhomogeneous pa- 
rameter ocj. For simplicity, we start by considering the periodic boundary conditions. 
The Hilbert space of the system is 

H= ®V;= (C 2 )^ N 

7=1 

In the QISM CHUB EH the so-called monodromy matrix 

T(A) = R m(A -«„)••■ Roi (A - «i) (III-l) 

is used to describe the system. For simplicity we have omitted the dependence on 
the quasi-classical parameter rj, the deformation parameter £ and the inhomogeneous 
parameters {ocj,j = 1, . . . ,N}. Notice that T(A) is a two-by-two matrix in the auxil- 
iary space Vq = C 2 , whose entries are operators acting in TL. Due to the Yang-Baxter 
equation (|II.12[) , it is straightforward to check that the monodromy matrix satisfies the 
RTT-relations IMEQUII 

R 12 (A - }i)T{\)T{ ¥ ) = T( M )T(A)R 12 (A - F ). (III.2) 
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The above equation is written in the tensor product of the auxiliary space Vq (g> Vq = 
C 2 (g) C 2 , using the standard notation of the QISM. 

The periodic boundary conditions and the RTT-relations (|III.2|l imply that the trans- 
fer matrix 

t(A) = tr T(A), (III.3) 

commute at different values of the spectral parameter, 

[t(A),*00]=0, (HI.4) 

here we have omitted the nonessential arguments. 

Due to the quasi-classical property dll.131 ) and the normalization of the R-matrix 
HII. 17^ . the quasi-classical expansion of the transfer matrix, for the special values of the 
spectral parameter, is given by 113711381 

Z k = t{A = * k ) = t + V H k + 0( V 2 ), (111.5) 

where H k are the corresponding Gaudin Hamiltonians, in the periodic case. The com- 
mutativity of the Gaudin Hamiltonians, 

[H k ,Hi}=0, (III.6) 

is ensured by the commutativity of the transfer matrix (|III.4|) and the fact the first term 
in the above expansion is the identity matrix, due to (III.13D and (|II.17[) . 

In order to construct integrable spin chains with non-periodic boundary condition, 
we use the Sklyanin formalism 11331 . The corresponding monodromy matrix 7~(A) con- 
sists of the two matrices T(A) flHLl[) and a reflection matrix K~(\) dll.301 ), 

T(A) = T(A)K-(A)T- 1 (-A) 

= R 0N (A - a N ) ■ ■ ■ R i(A - ai )K (\)R w (A + «i) • • • R N o(A + x N ), (III.7) 

where, for simplicity, we have suppressed the dependence on the other parameters. By 
construction, the exchange relations of the monodromy matrix 7"(A) in Vb <8> Vq are 

R 12 (A - ^)T(A)R 21 (A + fi)T(it) = T()*)Ri2(A + ^)T(A)R 21 (A - ft), (III.8) 

using the notation of |33). The open chain transfer matrix is given by the trace of T(A) 
over the auxiliary space Vo with an extra reflection matrix K + (\) 1331 , 

f(A)=tr (K + (A)T(A)). (III.9) 

The reflection matrix K + (A) dll.29l > is the corresponding solution of the dual reflection 
equation dll.281) . The commutativity of the transfer matrix for different values of the 
spectral parameter 

[t(\),t(ft)]=0, (111.10) 
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is guaranteed by the dual reflection equation dll.281 ) and the exchange relations ( 111181 ) of 
the monodromy matrix 7~(A). 

Analogously to the periodic case, the quasi-classical expansion of the transfer ma- 
trix, for special values of the spectral parameter, yields the corresponding Gaudin Hamil- 
tonians with boundary terms [40 1 



Z k = f (A = «*) = ! + rjH k + 0{rj 2 ) f (III.ll) 

where H^ are the Gaudin Hamiltonians. As in the periodic case, the commutativity of 
the Hamiltonians Hk is guaranteed by the equation dlll.lOb and the fact the first term 
in the above expansion is the identity matrix, due to the quasi-classical property of the 
R-matrix (1IL13I ) and the normalization conditions (|II.17[) , (1IL32I ) and dll.341 ). The Gaudin 
Hamiltonians will be calculated in the following section. 

IV Deformed Gaudin Hamiltonians with boundary terms 

In this section we calculate explicitly the Gaudin Hamiltonians. In the periodic case 
|37H38l, it is straightforward to calculate the first two terms in the expansion (1III.5I ). By 
setting 6 = — £?/ and A = a. k and using the normalization of the R-matrix (|II.17|) we 
obtain 

Zfc|r/=o = tr (RonK - «n) • • ■ Rok{*k -«*)••■ &0l(«fc - Oil)) | )? =o 
= tr (Ron (a* - «n) ■ • • Vok ■ ■ ■ Roi(ock ~ «i)) 1^=0 
= tr (V ok ) = 1. (IV.l) 

In order to derive the corresponding Gaudin Hamiltonians, in the second term we also 
use the quasi-classical property of the R-matrix dll.131) , 

dZk _ / dR i(cik — ocA \ 

H k = -3— |, 7 =o = }^ tr o I RoN{K k -K N ) V ok ■ ■ ■ R 01 {a k - oci) I \ v=0 

+ £ tr Uon(^ - «n) • • • Vok ■ ■ ■ dRol ^ k ~^ . . ■ R 0l ( Kk - ecj)) | 7=0 

= E r « (**-*!)• ( IV - 2 ) 

In terms of the local spin- 1 /^ operators Sy = -^- these Hamiltonians are 

H ^ 2 E [ ^zr- + e ( s * s / + - s / +s *)l . ( IV - 3 ) 



lf t k \ a k~oci 
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where ~ stands for equality up to a constant additive term. By implementation of 
the algebraic Bethe ansatz, it was shown in [43| that the spectra of these Hamiltonians 
coincide with the spectra of the s£(2) invariant model, the case when £ = 0. 

In the case of non-periodic boundary, we consider the left and right reflection matri- 
ces K~(A) and K + (A) given by flll.301 ) and ( III.29D , respectively together with 8 = — £)/. 
In order to obtain the expansion piLllj) in (IIII.9I) we specify A = a k and calculate the 
first term 

2*1)7=0 = fro (K + (ct k )T(<x. k )) l^o = tr (K^(oc k )R ON (oc k - cc N ) ■ ■ ■ V ok ■ ■ • Roi(«fc - «i) x 
xK (a k )R w (a. k + act) ■ ■ ■ R ko (2oc k ) ■ ■ ■ R N0 {a k + oc N )) \ n = 
= tr {K (-K k )T ok K (a k )) = K k (a k )K k (-a k ) = 1. (IV.4) 

In the last step above we have used the normalization flII.32D . The Gaudin Hamiltonians 
with boundary terms are related to the second term in the expansion (|III.11|) , 

H k = ^|,=o = tr {^^\^oV ok K (^ +tr fa(-*)P^ fa)***^^ 
+ J> (x o -(-^ 8R °' ( ^-^\ = o P 0fc X -(. fc )) 

+ L tr o K o (-*k)V ok ^- -1^=0^0 («*) 

l<k V ar l 

+ E tr o U Q {-* k )V ok K {* k ) dRlo{ l k + ai) \^) . (IV.5) 

In the derivation above we have used the fact that the left reflection matrix K~ (A) does 
not depend on the quasi-classical parameter tj, (1II33I) . Finally, the Gaudin Hamiltonians 
can be expressed in a more concise form 

H k = T k (tx k ) +^r H fe - oci) + Y2 K k ( a kh k (oc k + oci)K k {-oc k ), (IV.6) 

where r;j(A) is the corresponding classical r-matrix pi.!4[) and 

T k (cc k ) =K k -(cc k ) ( ^^ | 7 =o + tr (M~«fc) V ok r k0 (2a k )) J . (IV.7) 

Notice that the second term on the righthand side of the (|IV.6b coincides with the 
Gaudin Hamiltonians PV.2|) and that the first and the third term are the boundary terms 
depending on the reflection matrices K~ and K + . Consequently, the Hamiltonians flIV.61) 
are rational functions of the boundary parameters £, <p, xp and the deformation parame- 
ter £. 
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V Conclusions 

We have derived the deformed s£(2) Gaudin Hamiltonians with boundary terms. Start- 
ing from the SL (2) -invariant Yang R-matrix deformed by the Jordanian twist element 
for the universal enveloping algebra U(s£(2)) and generic solutions of the associated 
reflection equation and the dual reflection equation, we have obtained the deformed 
inhomogeneous spin-% XXX chain. We have shown that the quasi-classical expansion 
of the transfer matrix of the chain yields the corresponding Gaudin Hamiltonians with 
boundary terms. 
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